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Abstract 



Starting from the Maxwell's equations and without resort to the paraxial approxima- 
tion, we derive equations describing stationary (l+l)-dimensional beams propagating at an 
arbitrary direction in an optical crystal with cubic symmetry and purely quadratic (x^^^) 
nonlinearity. The equations are derived separately for beams with the TE and TM polar- 
izations. In both cases, they contain x^'^^ and cubic (x^^^) nonlinear terms, the latter ones 
generated via the cascading mechanism. The final TE equations and soliton solutions to 
them are quite similar to those in previously known models with mixed X^'^^'X^^^ nonlin- 
earities. On the contrary to this, the TM model is very different from previously known 
ones. It consists of four first-order equations for transverse and longitudinal components of 
the electric field at the fundamental and second harmonics. Fundamental-soliton solutions 
of the TM model are also drastically different from the usual x^"^^ solitons, in terms of the 
parity of their components. In particular, the transverse and longitudinal components of the 
electric field at the fundamental harmonic in the fundamental TM solitons are described, 
respectively, by odd and single-humped even functions of the transverse coordinate. Ampli- 
tudes of the longitudinal and transverse fields become comparable for very narrow solitons, 
whose width is commensurate to the carrier wavelength. 
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1 Introduction 



Optical crystals with the cubic symmetry are the basis of contemporary photonics. In 
particular, spatial optical solitons generated by the x^^"* (quadratic) nonlinearity, see reviews 
||T], 1^, are possible in these media. In the simplest case, these are two-component pulses 
composed of fundamental- frequency (FF) and second-harmonic (SH) waves. 

Usually, solitons in noncentrosymmetric crystals are considered when directions of the 
light propagation and polarization coincide with the crystal axes. However, in cubic crystals 
it is possible to investigate solitons carried by a beam which is propagating in arbitrary 
direction, the situation being similar to that known in electrooptics 0. In the present 
work, our objective is to study spatial solitons of both TE and TM (transverse-electric 
and transverse-magnetic) types with the arbitrary propagation direction of the carrier wave. 
First, we derive a corresponding general model (actually, we obtain a system of ODEs to 
describe a stationary shape of the (l+l)-dimensional solitons). The model is derived directly 
from the Maxwell's equations, without using the paraxial approximation, so that the final 
equations apply as well to very narrow (subwavelength) spatial solitons. 

In the case of the TE polarization, which is considered in section 2, the obtained model 
differs from the well-known simplest model giving rise to the x^^'' solitons by a single addi- 
tional cubic (x^^"*) term in the equation for the transverse FF component, which is generated 
by the underlying quadratic nonlinearities via the cascading mechanism (see a review |Q). As 
solitons in this TE model are very close to those in the known x'-^-'-x*'^'' models, we consider 
them, also in section 2, only briefly. 

Unlike the TE case, the final model for the TM solitons, which is derived in section 
3, turns out to be essentially novel, consisting of four first-order equations. Two equations 
for transverse FF and SH electric fields contain only quadratic nonlinear terms, while two 
other equations, for longitudinal electric fields, contain cascading-generated cubic terms too. 
Then, in section 4 we find, in a numerical form, stationary solutions for both fundamental 
and higher-order solitons in the TM model. The fundamental TM solitons turn out to 
be drastically different, as concerns the parity of their components, from solitons in the 
traditional x*^^'* models: their longitudinal FF electric field, which has no counterpart among 
components of the usual x^^'' solitons, is described by an even single-humped function, while 
the transverse TE electric field, which is an analog of the FF field in the usual solitons, is 
an odd function of the transverse coordinate x. Moreover, the transverse SH electric field, 
although being an even function, is not monotonically decaying in the interval < x < oo, 
but rather changes its sign at a finite x. 

The x*^^^ contribution to the polarization induced by the electric field with real vectorial 
components En in noncentrosymmetric dielectric crystals can be represented in the form 
P*^^) = dimnEmEn (summatiou over repeated subscripts is assumed), where dimn is a nonlinear 
susceptibility tensor 0. Obviously, it is symmetric with respect to the permutation of 
the subscripts m and n. For cubic noncentrosymmetric crystals in the crystollagraphic 
coordinates x.y^z, nonvanishing components of the tensor are di22, = (^213 = '^312 = 0- 

We will consider light propagation in an arbitrary direction, relative to the axes x,y,z, in 
a new coordinate frame (x, y, z), the 2;-axis being aligned with the propagation direction. The 
new frame may be generated by consecutive rotations of the initial one by angles (f, and 7 
around the z, x, and z axes, respectively. The matrix (a^n) of the orthogonal transformation 
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to the new frame, Xk = akiXi, where Xk and Xk stand for the sets (x, y, z) and (x, y, z), is 

cos (p cos 7 — cos {} sin ip sin 7 sin ip cos 7 + cos -d cos 99 sin 7 
— cos (y9 sin 7 — cos "i? sin ip cos 7 — sin (y9 sin 7 + cos cos ip cos 7 
— sin sin sin cos 99 

Components of the dielectric displacement vector in the new frame are 

Dk = sEk + 2dkmnEmEn, (2) 

where e is the dielectric constant, and 

dkmn = O'klO'mpO'nqdlpq = 2d{aklOm20n3 + Ok20mlOn3 + Ok30mlOn2) 

is the x*"^"* susceptibility tensor in the new frame. In the general case, all the 27 tensor 
components are different from zero, but they depend only on one constant d and three 
angles ip, and 7. 

Below, the full electric and magnetic fields, including both FF and SH components, will 
be sought for in the form 

Ek = afcCos(2<l>) +6fcsin(2<l>) +ylfcCos$ + -Bfc sin$, (3) 
Hk = PfcCos(2<l>) +gfcsin(2$) +PfcCos<l> + (5fcSin<l>. (4) 

Here, the FF phase is ^ = (3z — cut, u is the frequency, t is time, and {3 is the propagation 
constant, the small and capital letters standing for the SH and FF amplitudes, respectively. 
In the general case, with the carrier wave propagating in an arbitrary direction, solutions 
are of the mixed TE-TM type, i.e., they have all the electric and magnetic field components. 
However, pure TE and TM modes will be found below for certain angles ip, d and 7, including 
nontrivial cases when the propagation and polarization directions are not aligned with the 
crystallographic axes {x,y,z). 

To conclude the introduction, it is relevant to mention several other models of the x^^"* 
type which take into regard the polarization of the fields. They describe bimodal fields in 
an isotropic medium or the so-called type-II interactions between the FF and SH waves, 
with a single (see, e.g., Refs. 0) or double parametric resonance between them. However, 
none of these previously studied models included longitudinal components of the fields, nor 
magnetic field was explicitly introduced in them. 

2 The TE case 

First, we consider a simpler case with the FF wave of the TE type, i.e., when the electric 
FF field has only one component A2, which is a function of x, but directed along the y axis. 
Then, the only nonzero FF components are A2 = A, Pi = P, and Q3 = Q, the latter one 
being longitudinal magnetic field. Substituting Ek and Hk from Eqs. @ and (|) into the 
Maxwell's equations, it is possible to eliminate the magnetic-field components and obtain a 
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set of differential equations for tlie electric components only: 
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(5) 
(6) 

where the prime stands for d/ dx, jSo = uj/c is the FF propagation constant in vacuum, c being 
the light velocity in vacuum, /3f 2 = ^1,2^0, and £1,2 are FF and SH dielectric permeabilities. 
From the condition that the fields must vanish at x = ±00, it immediately follows from the 
last equation in the system that 62 = 0, and it can be shown too that gi = ^3 = 0. 

Note that, while /3 is the actual propagation constant for the FF component of the 
soliton-carrying beam, f3i and 2/^2 are the wave numbers of linear plane waves propagating 
in the medium at the frequencies u and 2uj. In fact, (3, which is different from f3i because of 
the nonlinearity and finite transverse size of the soliton beam, is a parameter of the soliton 
family at fixed values of the carrier-wave propagation constants /?i^2- 

In addition to the differential equations (|^), it also follows from the Maxwell's equations 
that the SH electric filed must satisfy a set of linear algebraic relations, dumO'm = '^32m«m = 
di2mbm = d22mbm = d^2mbm = 0. We do uot display details of analysis which shows that 
Eqs. (|^) are compatible with the extra linear relations at some special values of the angles 
(yj, "(9,7). After lengthy transformations, in all such cases Eqs. (H) can be reduced to a simple 
system of two equations. 



drf 
d^u 
drf 



R-uR + 6-R\ 

13^ - /j| 2 



(7) 
(8) 



where R and u are certain combinations of the FF and SH electric fields, r] = \J (3"^ — Plx, 
and 6 is some other constant. 

In particular, a nontrivial compatibility case can be obtained by choosing cos = cos 7 = 
cos(2<y5) = 0. One then has di22 = ±c?, "2 = "3 = &i = 0, and 



R = V2D {P/P2) A, u = D{ai + ^7^2), h = - 



du 



(9) 



with 6 = {P^- pf) 74/32, D = 2rfi22wVc. 

A special case when all the equations are compatible and S = (no effective cubic 
nonlinearity) corresponds (for instance) to cos(2'(9) = simp = sin 7 = 0, which again yields 
c?i22 = ±rf. In this case, Eqs. (0) and (H) are equivalent to the simplest version of the 
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stationary x^'^^ model, hence they have well-known soliton solutions |@|. If 5 7^ 0, the system 
differs from the simplest one by the extra x^^^ term in Eq. (^. 

Soliton-generating models combining x*-^-* and x^^^ nonlinearities are well known too, see, 
e.g., Refs. 0. Equations (|^ and (^) are somewhat different from them, as they contain the 
cubic term solely in the FF equation. In models of this type, soliton solutions can be easily 
found 0. A typical example of a TE fundamental (single- humped) soliton obtained from 
Eqs. (D, (D and (g) is displayed in Fig. 1. 

Cases when all the equations are compatible and, eventually, amount to the simplest 
system of Eqs. (|^) and (|) are possible not only for discrete sets of values of the three angles, 
as in the two examples shown above, but also for certain one-parameter families, when two 
equations are imposed on the three angles. Examples are 

sin (p = 2 sin 7, sin vcos 7 = 1; (10) 

cos^ (y9 = 2 sin^ 7, cos^9 = tan(y9 ■ tan7. (11) 

These examples do not exhaust all the compatibility cases. A complete analysis of the 
compatibility conditions for the TE type of the fields in the cubic crystal proves to be very 
tedious. 



3 The TM case 

We now proceed to a more interesting TM model, when nonzero FF amplitudes are Ai = 
A, A^, and P2 = P, see Eqs. (0) and (^. For this model, stationary equations for the 
magnetic- and electric-field amplitudes can be derived in a general form, which, however, 
turns out to be very cumbersome. Therefore, we here concentrate on a particular (but quite 
nontrivial) case, 

sin2y) = sin 7 = cos2'(9 = 0. (12) 
In the present case, equations for nonvanishing fields can be cast into the following form 

(3P = -(3oieiA + 2dA^h), 

P' = (3o[eiA3 + 2d{Ah - A;a,)l 

PoP = -f3A-A',, 

(3p2 = -Po (^e2ai - dAfj , 

p'2 = 2(3o{e2b3 + 2dAA3), 

2P0P2 = -2/3a^~b',. (13) 

It is convenient to define renormalized transverse {Ez,ez) and longitudinal {E,e) FF and SH 
electric fields, 

A, . -E., 'Ma^E, (14) 



/32 - (31(3 



2d(3l , _ 2d(3l _ 

^ " -ai = e. (15) 



/?2 - PIP P^ 
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Finally, the magnetic fields P and p2 can be eliminated from Eqs. ([T3|) , which yields a final 
fourth-order system for the electric fields (where, again, i] = \J P"^ — Plx), 



E - E,e,, (16) 



dE^ 
drj 

This model with mixed x^^"* ^i^iid x*^^) nonlinearities appears to be novel, as compared to 
various known models of the X^'^^'X^^^ typ^ 0- 

It is relevant to stress that, since both the TE and TM models, composed of Eqs. (|^), 
(§) and (|16|), ([T7|), respectively, have been derived directly from the Maxwell's equations, 
without resorting to the paraxial approximation 0, the applicability of these models is not 
limited to the case of broad solitons, whose width is much larger than the carrier wavelength. 
In other words, the final equations derived above apply as well to very narrow subwavelength 
solitons, which have recently attracted considerable attention in spatial and temporal models 
with the Kerr nonlinearity [1^], but were not thus far considered in x^^^ systems. 



4 TM solitons 

As the system of equations (|16D and (|l^ is essentially different from the standard models, 
it is necessary to investigate its soliton solutions, which we have done by means of the usual 
shooting technique. Fundamental and higher-order solitons were found in a broad parametric 
region. 

The first noteworthy result is that fundamental solitons are characterized by a single- 
humped profile of the longitudinal electric FF field, corresponding to an even single-humped 
function Ez{x), while the transverse FF field is described by an odd function E{z). In other 
words, only the distribution of the squared longitudinal electric FF field, -E^(x), is single- 
humped in the fundamental TM soliton, while the distribution of the squared transverse FF 
field, E^{x), is double-humped, with E'^{0) = 0. As for the SH components, it is evident 
from Eqs. ([T7| ) that the corresponding transverse and longitudinal fields e and must be, 
respectively, even and odd functions of x, whatever parity of the FF fields E and E^. This 
feature is confirmed by Fig. 2; however, a nontrivial property of the fundamental soliton is 
that the even function e{x) is not monotonic in the interval < x < oo, giving rise to three 
humps and two zeros in the distribution of the corresponding squared field e^(x). 

Note that the opposite parities of the FF transverse and longitudinal fields E{x) and 



Ez{x) is an obvious consequence of Eqs. (pIBI). Evidently, besides the structure seen in 



Fig. 2, with the odd transverse and even longitudinal FF fields, the latter condition is also 
compatible with the opposite case, when E{x) and Ez{x) would be, respectively, even and 
odd. However, numerical solutions have never turned up TM solitons of that type. 
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Lastly, Fig. 2 clearly suggests a conclusion which is also supported by numerical solutions 
obtained at many other values of the parameters: in terms of the renormalized variables 
defined by Eqs. (0) and (H), the maximum values of the longitudinal and transverse fields 
are nearly equal. To compare the maximum values of the physical fields, one should undo 
the rescalings ([l^) and (1151) . An obvious result is that the amplitudes of the longitudinal 
FF and SH electric fields remain much smaller than the amplitudes of the transverse fields 



if 



f3i <^ (3, which is the case for the usual solitons whose width is much larger than 
the carrier wavelength. For narrow solitons, whose width is comparable to or smaller than 
the wavelength, the physical amplitudes of the longitudinal and transverse fields become 
comparable too. 

Thus, the properties of the TM soliton are drastically different from, and essentially 
more nontrivial than those of the TE soliton. Indeed, the latter one has the transverse FF 
electric field described by an even function ^2(0;), which is monotonically decaying in the 
interval < a; < 00. Another strong difference is that, although the SH transverse electric 
field in both the TM and TE solitons is described by an even function, only in the latter 
case this function is monotonic at < x < 00. 

It is also natural to compare the TE and TM solitons found above with solitons in the 
usual x^"^^ models, which, as a matter of fact, are also of the TE type. As well as the TE 
soliton shown in Fig. 1, fundamental solitons in the usual models are characterized by even 
distributions of both FF and SH electric fields, with the single hump at the center of the 
soliton. Thus, the fundamental TM solitons found in this work are drastically different from 
the traditional x*^^"* solitons. 

The numerical integration of Eqs. (0) and (|l^) readily generates, alongside the fun- 
damental TM solitons, higher-order ones. A typical example of a second-order TM soliton, 
in which each component has an additional extremum inside the interval < x < 00, is 
displayed in Fig. 3. It is quite easy to find solitons of still higher orders too. 

It is relevant to mention that higher-order solitons are well known in usual x^"^^ models, 
both in the simplest two-component systems (see, e.g., Refs. |jll|) and in more sophisticated 
three- [|12| and four-component models. In most cases, higher-order solitons are subject 
to strong dynamical instability. Nevertheless, they were found to be (numerically) stable in 
three- and four-component models combining the x^"^^ nonlinearity and resonant reflections 
on a Bragg grating |jl3| . Direct analysis of the soliton stability in the present system is quite 
difficult and is beyond the scope of this work, as it is necessary to simulate the full system 
of the Maxwell's equations. Nevertheless, it seems very plausible that both the TE and TM 
fundamental solitons may be stable, while the higher-order ones are unstable. 



5 Conclusion 

We have derived, starting from the Maxwell's equations and without resort to the paraxial 
approximation, equations describing stationary (l-l-l)-dimensional beams propagating at an 
arbitrary direction in an optical crystal with cubic symmetry and purely quadratic x^^"* 
nonlinearity. The equations were derived separately for beams of the TE and TM types. 
In both cases, they contain x^^^ and x*"^"* nonlinear terms, the latter ones being generated 
via the cascading mechanism. The final TE equations and soliton solutions to them are 
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quite similar to those in previously known models with mixed x^'^^-X^^^ nonlinearities. On 
the contrary to this, the TM model turns out to be very different from the previously known 
ones. It consists of four first-order equations for transverse and longitudinal components of 
the electric field at the fundamental and second frequencies. Fundamental-soliton solutions 
to the TM model are also drastically different from the usual x^^-* solitons, in terms of the 
parity of their components. Most noteworthy, the transverse and longitudinal components 
of the electric field at the fundamental frequency in these solitons have, respectively, odd 
and single-humped even shapes. Besides that, it was also concluded that maximum values 
of the longitudinal and transverse fields become comparable for very narrow solitons, whose 
width is commensurate to or smaller than the carrier wavelength. 



8 



References 

[1] M. Segev, Opt. Quant. Electron. 30, 503 (1998); Y.S. Kivshar, Opt. Quant. Electron. 
30, 571 (1998). 

[2] C. Etrich, F. Lederer, B. Malomed, T. Peschel, and U. Peschel, Progr. Opt., vol. 41 (in 
press) . 

[3] B.V. Gisin, Kristallografiya 16, 151 (1971) [Sov. Phys. Crystallogr. 16, 118 (1971)]; 
D.F. Nelson, J. Opt. Soc. Am. 65, 1144 (1975). 

[4] G.l. Stcgcman, D.J. Hagan, and L. Torncr, Opt. Quant. Elcctr. 28, 1691 (1996). 

[5] N. Blombergen, Nonlinear Optics (Benjamin: New York, 1965); A. Yariv and P. Yeh, 
Optical waves in crystals (John Willey & Sons: New York, 1984). 

[6] A.D. Boardman, P. Bontemps, and K. Xie, Opt. and Quant. Electron. 30, 891 (1998). 

[7] H.T. Tran, Opt. Comm. 118, 581 (1995); D. Anderson, M. Lisak, and B.A. Malomed, 
Opt. Comm. 126, 251 (1996); A.V. Buryak, Y.S. Kivshar, and S. Trillo, J. Opt. Soc. 
Am. B 14, 3110 (1997); U. Peschel et al., Phys. Rev. E 55, 7704 (1997). 

[8] S. Trillo and G. Assanto, Opt. Lett. 19, 1825 (1994); I. Towers, A.V. Buryak, R.A. 
Sammut, and B.A. Malomed, Opt. Lett. 24, 1738 (1999). 

[9] S. Trillo, A.V. Buryak, and Y.S. Kivshar, Opt. Comm. 122, 200 (1996); O. Bang, Y.S. 
Kivshar, and A.V. Buryak, Opt. Lett. 22, 1680 (1997); J. Yang, B.A. Malomed, and 
D.J. Kaup, Phys. Rev. Lett. 83, 1958 (1999) . 

[10] A.W. Snyder, D.J. MitcheU, and Y. Chen, Opt. Lett. 19, 524 (1994); C.V. Hile and 
W.L. Kath, J. Opt. Soc. Am. B 13, 1135 (1996); E.Granot et al.. Opt. Lett. 22, 1290 
(1997); Opt. Comm. 166, 121 (1999); ibid. 178, 431 (2000); C. Chen and S. Chi, Opt. 
Commun. 157, 170 (1998); S. Blair, Chaos 10, 570 (2000). 

[11] A.V. Buryak and Y.S. Kivshar, Phys. Lett. A 197, 407 (1995); H. He, M.J. Werner, 
and P.D. Drummond, Phys. Rev. E 54, 896 (1996). 

[12] W.C.K. Mak, B.A. Malomed, and P.L. Chu, Phys. Rev. E 58, 6708 (1998); A.R. Champ- 
neys and B.A. Malomed, Phys. Rev. E 61, 886 (2000). 

[13] T. Peschel et al., Phys. Rev. E 55, 4730 (1997). 



9 



Figure Captions 

Fig. 1. A typical example of a bright spatial soliton with the TE polarization at values of 
the parameters — /?|) / — (3f) = 0.8 and S = 0.1. Shown in this figure are the fields 
R and u, obtained from Eqs. (7) and (8), vs. the normalized transverse coordinate rj. 

Fig. 2. A typical example of a fundamental bright spatial sohton with the TM polariza- 
tion. Here and in the subsequent figure, the continuous and dashed lines show, respectively, 
the transverse and longitudinal normalized fields: (a) the FF fields E and E^, and (b) the 
SH fields e and 63. The parameters are = 0.60, (/52//5)^ = 0.80. 

Fig. 3. A typical example of a second-order bright spatial soliton with the TM polar- 
ization. The parameters are = {^2/^)'^ — 0.35. 
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